Measurement of Muon Neutrino Quasi-Elastic Scattering on Carbon 



A. A. Aguilar-Arevalo 5 , A. 0. Bazarko 12 , S. J. Brice 7 , B. C. Brown 7 , L. Bugcl 5 , J. Cao 11 , L. Coney 5 , 
J. M. Conrad 5 , D. C. Cox 8 , A. Curioni 16 , Z. Djurcic 5 , D. A. Finlcy 7 , B. T. Fleming 16 , R. Ford 7 , F. G. Garcia 7 , 
G. T. Garvey 9 , C. Green 7 ' 9 , J. A. Green 8 - 9 , T. L. Hart 4 , E. Hawker 15 , R. Imlay 10 , R. A. Johnson 3 , P. Kasper 7 , 
T. Katori 8 , T. Kobilarcik 7 , I. Kourbanis 7 , S. Koutsoliotas 2 , E. M. Laird 12 , J. M. Link 14 , Y. Liu 11 , Y. Liu 1 , 
W. C. Louis 9 , K. B. M. Mahn 5 , W. Marsh 7 , P. S. Martin 7 , G. McGregor 9 , W. Mctcalf 10 , P. D. Meyers 12 , 
F. Mills 7 , G. B. Mills 9 , J. Monroe 5 , C. D. Moore 7 , R. H. Nelson 4 , P. Nienabcr 13 , S. Oucdraogo 10 , 
R. B. Patterson 12 , D. Pcrcvalov 1 , C. C. Polly 8 , E. Prebys 7 , J. L. Raaf 3 , H. Ray 9 , B. P. Roe 11 , A. D. Russell 7 , 
V. Sandberg 9 , R. Schirato 9 , D. Schmitz 5 , M. H. Shacvitz 5 , F. C. Shoemaker 12 , D. Smith 6 , M. Sorel 5 , 
P. Spcntzouris 7 , I. Stancu 1 , R. J. Stefanski 7 , M. Sung 10 , H. A. Tanaka 12 , R. Tayloe 8 , M. Tzanov 4 , 
R. Van dc Water 9 , M. O. Wascko 10 , D. H. White 9 , M. J. Wilking 4 , H. J. Yang 11 , G. P. Zcllcr 5 , E. D. Zimmerman 4 

(The MiniBooNE Collaboration) 
1 University of Alabama; Tuscaloosa, AL 35487 
2 Bucknell University; Lewisburg, PA 17837 
3 University of Cincinnati; Cincinnati, OH 45221 
4 University of Colorado; Boulder, CO 80309 
5 Columbia University; New York, NY 10027 
6 Embry-Riddle Aeronautical University; Prescott, AZ 86301 
7 Fermi National Accelerator Laboratory; Batavia, IL 60510 

^Indiana University; Bloomington, IN 47405 
9 Los Alamos National Laboratory; Los Alamos, NM 87545 
10 Louisiana State University; Baton Rouge, LA 70803 
11 University of Michigan; Ann Arbor, MI 48109 
12 Princeton University; Princeton, NJ 08544 
13 Saint Mary's University of Minnesota; Winona, MN 55987 
14 Virginia Polytechnic Institute & State University; Blacksburg, VA 24061 
Western Illinois University; Macomb, IL 61455 
16 Yale University; New Haven, CT 06520 

(Dated: February 1, 2008) 

The observation of neutrino oscillations is clear evidence for physics beyond the standard model. 
To make precise measurements of this phenomenon, neutrino oscillation experiments, including 
MiniBooNE, require an accurate description of neutrino charged current quasi-elastic (CCQE) cross 
sections to predict signal samples. Using a high-statistics sample of CCQE events, MiniBooNE 
finds that a simple Fermi gas model, with appropriate adjustments, accurately characterizes the 
CCQE events observed in a carbon-based detector. The extracted parameters include an effective 
axial mass, = 1.23 ± 0.20 GeV, that describes the four-momentum dependence of the axial- 

vector form factor of the nucleon; and a Pauli-suppression parameter, k = 1.019 ± 0.011. Such 
a modified Fermi gas model may also be used by future accelerator-based experiments measuring 
neutrino oscillations on nuclear targets. 



The recent observation of neutrino oscillations is strong 
evidence for massive neutrinos and, therefore, for physics 
beyond the standard model. Accelerator-based experi- 
ments searching for neutrino oscillations, such as Mini- 
BooNE [lj and K2K Q , use charged current quasi-elastic 
(CCQE) interactions to search for the appearance of elec- 
tron neutrinos (y e n — > e~p) in beams of muon neutrinos. 
The muon neutrino CCQE interaction (v^n — > pL~p) thus 
provides a calibration for the neutrino beam and for the 
interaction cross section. In addition, such events domi- 
nate at energies between 200-2000 MeV where the oscilla- 
tion searches are conducted. To ensure high event yields, 
these experiments use nuclear media (carbon or water) 
as the neutrino target; therefore, it is crucial to employ 
an accurate model of the CCQE interaction on nuclei. 



In this Letter, we describe the model improvements de- 
veloped for the recent oscillation search from the Mini- 
BooNE experiment [lj|. The modified model describes 
this reaction remarkably well and should be relevant for 
future accelerator-based neutrino oscillation searches. 

To model the scattering from nucleons confined in nu- 
clei, most neutrino oscillation experiments employ an 
event generator based on the relativistic Fermi gas (RFG) 
model Q ■ Such models assume a flat nucleon momentum 
distribution up to some Fermi momentum (j>f), assign a 
single value for the nucleon binding energy (Eb) to ac- 
count for the initial and final state total energies, and 
utilize standard nucleon vector and axial-vector on-shell 
form factors. Many of these model parameters may be 
inferred from existing data; for example, pf, Eb, and the 
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vector form factors can be determined from elastic elec- 
tron scattering data [1, Q . Despite providing these con- 
straints, electron data yield limited information on the 
axial-vector form factor of the nucleoli and the CCQE 
cross section at very low four-momentum transfer (Q 2 ). 
Present knowledge of the axial-vector form factor has 
been informed largely by past neutrino experiments, but 
these suffer from low statistics and were performed using 
predominantly deuterium targets [f|. Since these early 
measurements, neutrino experiments have encountered 
difficulties describing their data at low Q 2 , where nuclear 
effects arc largest, and have often measured axial- vector 
form factor parameters above some minimum Q 2 value. 

The MiniBooNE experiment has collected the largest 
sample of low energy muon neutrino CCQE events to 
date. We describe here the use of such events in tuning 
the REG model to better describe quasi-elastic scattering 
on nuclear targets. The analysis fits the reconstructed Q 2 
distribution of the MiniBooNE CCQE data in the region 
< Q 2 < 1 GeV 2 to a simple RFC model Q with two 
adjustable parameters: the axial mass, Ma, appearing in 
the axial- vector form factor; and k, a parameter that ad- 
justs the level of Pauli-blocking at low values of Q 2 . The 
best-fit model results in a good description of the data 
across the full kinematic phase space including the low- 
Q 2 region. This technique is crucial to the MiniBooNE 
oscillation search [I[ as it is used to predict the v e CCQE 
oscillation events based on the constraints provided by 
the high-statistics MiniBooNE CCQE sample. 

The Fcrmilab Booster neutrino beam, optimized for 
the MiniBooNE oscillation search, is particularly suited 
for investigation of low energy neutrino interactions. The 
Fermilab Booster provides 8.89 GeV/c protons which col- 
lide with a 71 cm long beryllium target inside a mag- 
netic horn. The horn focuses positively charged pions 
and kaons produced in these collisions, which can subse- 
quently decay in a 50 m long decay region, yielding an in- 
tense flux of muon neutrinos. A GEANT4-based [7| beam 
simulation uses a parametrization Q of pion production 
cross sections based on recent measurements from the 
HARP Q and E910 experiments, along with a de- 
tailed model of the bcamlinc geometry to predict the neu- 
trino flux as a function of neutrino energy and flavor. The 
resulting flux of neutrinos at the MiniBooNE detector is 
predicted to be 93.8% (5.7%) (D^) with a mean en- 
ergy of ~ 700 MeV. Because 99% of the flux lies below 
2.5 GeV, the background from high multiplicity neutrino 
interactions is small. Approximately 40% of the total 
events at MiniBooNE are predicted to be CCQE, of 
which 96% result from pion decays in the beam. 

The MiniBooNE detector is a spherical tank of inner 
radius 610 cm filled with 800 tons of mineral oil (CH2), 
situated 541 meters downstream of the proton target. An 
optical barrier divides the detector into two regions, an 
inner volume with a radius of 575 cm and an outer vol- 
ume 35 cm thick. The inner region of the tank houses 



1280 inward-facing 8 inch photomultiplier tubes (PMTs), 
providing 10% photocathodc coverage. The outer region 
is lined with 240 pair-mounted PMTs which provide a 
veto for charged particles entering or leaving the tank. 
Muons produced in CCQE interactions emit primarily 
Cherenkov light with a small amount of scintillation light. 
A large number of muons stop and decay in the main de- 
tector volume. The muon kinetic energy resolution is 
7% at 300 MeV and the angular resolution is 5°. The 
response of the detector to muons is calibrated using a 
dedicated muon tagging system that independently mea- 
sures the muon energy for cosmic ray muons ranging up 
to 800 MeV. 

Neutrino interactions within the detector are simu- 
lated with the v3 NUANCE event generator [ll|. This 



program provides the framework for tuning the CCQE 
cross section parameters (described below) and predicts 
backgrounds to the sample, including neutrino induced 
single pion production events (CC l7r). Pion interac- 
tions in the nucleus and photon emission from nuclear 
de-excitation in nuance are tuned to reproduce Mini- 
BooNE and other [H data. A GEANT3-based El de- 



tector model (with GCALOR [14j hadronic interactions) 
simulates the detector response to particles produced in 
neutrino interactions. The simulation of light production 
and propagation in mineral oil has been tuned using ex- 
ternal small-sample measurements [ll| , muon decay elec- 
trons (also used to calibrate the energy scale) , and recoil 
nucleons from neutrino neutral current (NC) elastic scat- 
tering events. The predicted events arc additionally over- 
laid with events measured in a beam-off gate, in order to 
incorporate backgrounds from natural radioactivity and 
cosmic rays into the simulated data. 

Because of the low energy neutrino beam and Mini- 
BooNE detector capabilities, the identification of 
CCQE interactions relies solely on the detection of the 
primary muon and associated decay electron in these 
events: 
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This simple selection is highly effective for several rea- 
sons. First, the efficiency for detecting the decay of the 
produced in such events is high, 83%. The losses are 
due to muon capture on carbon (8% [I3|) and insuffi- 
cient decay time or energy of the decay electron (10%). 
Second, the CC l7r + contamination is significantly re- 
duced by requiring a single decay electron, since CC l7T + 
events typically yield two decay electrons, one each from 
the primary muon and the 7r + decay chains. The ex- 
ceptions are cases in which the primary \jT is captured 
or, more likely, the 7r + is either absorbed or undergoes 
a charge-changing interaction in the target nucleus or 
detector medium. Each of these processes is included 
in the detector simulation. Finally, by avoiding require- 
ments on the outgoing proton kinematics, the selection 
is inherently less dependent on nuclear models. 



Timing information from the PMTs allows the light 
produced by the initial neutrino interaction (first "sub- 
event") to be separated from light produced by the de- 
cay electron (second sub-event). The time and charge 
response of the PMTs is used to reconstruct the posi- 
tion, kinetic energy, and direction vector of the primary 
particle within each sub-event. Once separated into sub- 
events, we require that the first sub-event (the neutrino 
interaction) must occur in coincidence with a beam pulse, 
have a reconstructed position < 500 cm from the center 
of the detector, possess < 6 veto-PMT hits to ensure 
containment, and have > 200 main-PMT hits to avoid 
electrons from cosmic ray muon decays. The second sub- 
event (the \i~ decay electron) must have < 6 veto-PMT 
hits and < 200 main-PMT hits. Subsequent cuts specif- 
ically select CCQE events and discriminate against 
CC l7r + backgrounds. First, events must contain exactly 
two sub-events. Second, the distance between the elec- 
tron vertex and muon track endpoint must be less than 
100 cm, ensuring that the decay electron is associated 
with the muon track. 

A total of 193,709 events pass the MiniBooNE 
CCQE selection criteria from 5.58 x 10 20 protons on tar- 
get collected between August 2002 and December 2005. 
The cuts are estimated to be 35% efficient at selecting 
CCQE events in a 500 cm radius, with a CCQE purity of 
74%. The 35% efficiency is the product of a 50% prob- 
ability for containing events within the tank, the afore- 
mentioned 83% muon decay detection efficiency, and an 
85% efficiency for the electron vertex to muon endpoint 
requirement. 

The predicted backgrounds arc: 75% CC 1tt+, 15% CC 
1tt°, 4% NC 1tt ± , 3% CC multi-vr, 1% NC elastic, 1% 
CC ITT", 1% NC Itt , < 1% rj/p/K production, and < 1% 
deep inelastic scattering (DIS) and other events [ll|. In 
the analysis, cross section uncertainties of 25%, 40%, and 
25% are assumed on the Itt, mufti- 7r plus rj/p/K produc- 
tion, and DIS backgrounds, respectively. Because pions 
can be absorbed via final state interactions in the target 
nucleus, a large fraction of the background events look 
like CCQE events in the MiniBooNE detector. "CCQE- 
like" events, all events with a muon and no pions in the 
final state, are predicted to be 84% of the sample after 
cuts. 

The observables in the MiniBooNE CCQE sample 
are the muon kinetic energy T M , and the muon angle 
with respect to the neutrino beam direction 8^. The 
high-statistics MiniBooNE data sample allows us to ver- 
ify the simulation in two dimensions. Figure [T] shows the 
level of agreement between the shape of the data and 
simulation in the CCQE kinematic quantities before any 
CCQE cross section model adjustments. For this com- 
parison, the simulation assumes the RFC model as im- 
plemented in nuance @, EH, with E B = 34 MeV [H, 
Pf = 220 MeV/c updated non-dipole vector form 
factors [B| , and a non-zero pseudoscalar form factor 17 1 . 
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(a) E v =0.4GeV 

(b) E v =0.8GeV 

(c) E^l^GeV 

(d) Q 2 =0.2GeV 2 

(e) Q 2 =0.6GeV 2 

(f) Q 2 =1.0GeV 2 
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FIG. 1: Ratio of MiniBooNE CCQE data/simulation as a 
function of reconstructed muon angle and kinetic energy. The 
prediction is prior to any CCQE model adjustments and has 
been normalized to the data. The xVdof = 79.5/53. The 
ratio forms a 2D surface whose values are represented by the 
gray scale, shown on the right. If the simulation modeled the 
data perfectly, the ratio would be unity everywhere. Contours 
of constant E v and Q 2 are overlaid, and only bins with > 20 
events in the data are plotted. 



The axial-vector form factor is assumed to have a dipole 
form function of Q 2 with one adjustable parame- 
ter, Ma, the so-called "axial mass", Fa(Q 2 ) = 9a/ + 

q 2 /m\y. 

The simulation shown in Fig. [T] specifically assumes 
g A = 1.2671 [13 and M A = 1.03 GeV [H]. These model 
parameters are common defaults in most neutrino simu- 
lations. The figure shows that the disagreement between 
data and simulation follows lines of constant Q 2 and not 
E v . This supports the assumption that the data/model 
disagreement is not due to a mis-modeling of the incom- 
ing neutrino energy spectrum but an inaccuracy in the 
simulation of the CCQE process itself. We also explicitly 
assume no disappearance due to oscillations. 

Guided by indications that the data-model discrepancy 
is only a function of Q 2 , we have modified the existing 
CCQE model rather than introduce more drastic changes 
to the cross section calculation. This approach works well 
and requires adjustment of only two parameters: Ma and 
E\ Q . The parameter E\ Q effectively controls the effect of 
Pauli-blocking. It is the lower bound of integration over 
initial state nucleon energy and appears within the RFC 
model together with an upper bound E^: 



E hi = yjfi, + Ml E lo = ^p 2 F + M 2 - uj + E B , (1) 

where M n is the target neutron mass, M p is the outgo- 
ing proton mass, and u> is the energy transfer. In the 
RFC model, E^ is the energy of an initial nucleon on 
the Fermi surface and E\ is the lowest energy of an ini- 
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tial nucleoli that leads to a final nucleon just above the 
Fermi momentum (and thus obeying the exclusion prin- 
ciple in the final state). In practice, a simple scaling of 
E\ a was implemented in the MiniBooNE CCQE data fit 

via E\ = n(^Jp 2 F + M 2 — lu + E B ). The parameter k 

adds a degree of freedom to the RFG model which can 
describe the smaller cross section observed in the data at 
low momentum transfer and is likely compensating for 
the naive treatment of Pauli-blocking in the RFG model. 

The adjusted RFG model is then fit to the shape of 
the reconstructed Q 2 distribution in the MiniBooNE 
CCQE data: 

Q 2 = -q 2 = -ml + 2E V {E I1 - p„ cos 9„) > 0, (2) 

where is the muon mass, E^ (p^) is the reconstructed 
muon energy (momentum), and 8^ is the reconstructed 
muon scattering angle. The reconstructed neutrino en- 
ergy E v is formed assuming the target nucleon is at rest 
inside the nucleus: 

_ 2(M» - E B )E fl - (E% - 2M n E B + m* + AM 2 ) 

2 [(M„ -E B )-E ll + Pfl cos e„] ^ 

where AM 2 = M 2 — M 2 and E B > 0. A small correction 
is applied to E v in both data and simulation to account 
for the biasing effects of Fermi-smearing. This procedure, 
while yielding a more accurate E v estimate, has a neg- 
ligible impact on the Q 2 fit to MiniBooNE CCQE data. 
These expressions, with reconstructed muon kinematics, 
yield an E v resolution of 11% and a Q 2 resolution of 21% 
for CCQE events. 

The model parameters Ma and k are obtained from a 
least-squares fit to the measured data in 32 bins of recon- 
structed Q 2 from to 1 GcV 2 . All other parameters of 
the model are held fixed to the values listed previously, 
and a complete set of correlations between systematic un- 
certainties is considered. The total prediction is normal- 
ized to the data for each set of parameter values. Thus, 
the procedure is sensitive only to the shape of the Q 2 dis- 
tribution, and any changes in the total cross section due 
to parameter variation do not impact the quality of fit. 
The Q 2 distributions of data and simulation before and 
after the fitting procedure are shown in Figure [5] The 
X 2 /dof of the fit is 32.8/30 and the parameters extracted 
from the MiniBooNE CCQE data are: 

Mf = 1.23 ±0.20 GeV; (4) 
k = 1.019 ±0.011. (5) 

While normalization is not explicitly used in the fit, the 
new model parameters increase the predicted rate of 
CCQE events at MiniBooNE by 5.6%. The ratio of de- 
tected events to predicted, with the best-fit CCQE model 
parameters, is 1.21 ± 0.24. 

In general, varying Ma allows us to reproduce the high 
Q 2 behavior of the observed data events. A fit for Ma 
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FIG. 2: Reconstructed Q 2 for CCQE events including sys- 
tematic errors. The simulation, before (dashed) and after 
(solid) the fit, is normalized to data. The dotted (dot-dash) 
curve shows backgrounds that are not CCQE (not "CCQE- 
like"). The inset shows the la CL contour for the best-fit 
parameters (star), along with the starting values (circle), and 
fit results after varying the background shape (triangle). 

above Q 2 > 0.25 GcV 2 yields consistent results, M A S = 
1.25±0.12 GeV. However, fits varying only Ma across the 
entire Q 2 range leave considerable disagreement at low 
Q 2 (x 2 /dof = 48.8/31). The Pauli-blocking parameter k 
is instrumental here, enabling this model to match the 
behavior of the data down to Q 2 = (Figure [2]). 

Figure [3] shows the agreement between data and sim- 
ulation after incorporation of the Ma and k values from 
the Q 2 fit to MiniBooNE data. Comparing to Figure [TJ 
the improvement is substantial and the data are well- 
described throughout the kinematic phase space. 

Table [T] shows the contributions to the systematic un- 
certainties on Ma and k. The detector model uncertain- 
ties dominate the error in Ma due to their impact on 
the energy and angular reconstruction of CCQE events 
in the MiniBooNE detector. The dominant error on k 
is the uncertainty in the Q 2 shape of background events. 
This error (not included in the contour of Figure ^ is 
evaluated in a separate fit, where MiniBooNE CC l7r + 
data are used to set the background instead of the event 
generator prediction, and then added in quadrature. 

The result reported here, Mf = 1.23 ± 0.20 GeV, 
is consistent with a recent K2K measurement on a wa- 
ter target, M A = 1.20 ± 0.12 GeV Both values 
are consistent with but higher than the historical value, 
Ma = 1.026±0.021 GeV, set largely by deuterium-based 
bubble chamber experiments [l9j. The Ma value re- 
ported here should be considered an "effective param- 
eter" in the sense that it may be incorporating nuclear 
effects not otherwise included in the RFG model. In par- 
ticular, it may be that a more proper treatment of the 
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eter, k = 1.019 ±0.011, achieving substantially improved 
agreement with the observed kinematic distributions in 
1-15 this data set. Incorporation of both fit parameters al- 
lows, for the first time, a description of neutrino CCQE 
scattering on a nuclear target down to Q 2 = GeV 2 . 
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FIG. 3: Ratio of data/simulation as a function of muon ki- 
netic energy and angle after CCQE model adjustments. The 
simulation has been normalized to the data. The x 2 /dof = 
45.1/53. Compare to Figure [T] 



TABLE I: Uncertainties in M\ and k from the fit to Mini- 
BooNE Vy, CCQE data. The total error is not a simple 
quadrature sum due to correlations between the Q 2 bins cre- 
ated by the systematic uncertainties. 



error source 


SMf 


5k 


data statistics 


0.03 


0.003 


neutrino flux 


0.04 


0.003 


neutrino cross sections 


0.06 


0.004 


detector model 


0.10 


0.003 


CC 7r + background shape 


0.02 


0.007 


total error 


0.20 


0.011 



nucleon momentum distribution in the REG would yield 
an Ma value in closer agreement to that measured on 
deuterium. Future efforts will therefore explore how the 
value of Ma extracted from the MiniBooNE data is al- 
tered upon replacement of the RFG model with more 
advanced nuclear models 
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In summary, modern quasi-clastic scattering data on 
nuclear targets are revealing the inadequacies of present 
neutrino cross section simulations. Taking advantage of 
the high-statistics MiniBooNE CCQE data, we have 
extracted values of an effective axial mass parameter, 
Mf = 1.23 ± 0.20 GeV, and a Pauli-blocking param- 
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